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When a three-dimensional (3D) ferromagnetic topological insulator thin film is magnetized out-of- 
plane, conduction ideally occurs through dissipation-less, one-dimensional (ID) chiral states that are 
characterized by a quantized, zero-field Hall conductance. The recent realization of this phenomenon 
- the quantum anomalous Hall effect - provides a conceptually new platform for studies of edge- 
state transport, distinct from the more extensively studied integer and fractional quantum Hall 
effects that arise from Landau level formation. An important question arises in this context: how 
do these ID edge states evolve as the magnetization is changed from out-of-plane to in-plane? We 
examine this question by studying the field-tilt driven crossover from predominantly edge state 
transport to diffusive transport in Cr x (Bi,Sb) 2 - x Te 3 thin films, as the system transitions from a 
quantum anomalous Hall insulator to a gapless, ferromagnetic topological insulator. The crossover 
manifests itself in a giant, electrically tunable anisotropic magnetoresistance that we explain using 
the Landauer-Biittiker formalism. Our methodology provides a useful means of quantifying edge 
state contributions to transport in temperature and chemical potential regimes far from perfect 
quantization. 


In the absence of magnetic doping, surface electrons in 
3D topological insulators such as (Bi,Sb) 2 Te 3 are helical 
Dirac fermions protected from back scattering by time re¬ 
versal symmetry mm. When such a 3D topological insu¬ 
lator is made ferromagnetic by doping with magnetic ions 
pang or by integration into magnetic heterostructures [ 101 - 
EO, time-reversal symmetry is broken. If the magnetiza¬ 
tion M is perpendicular to the surface, the dispersion is 
gapped at the Dirac point [T31 Hi] , while M parallel to 
the surface is expected to keep the Dirac point intact 
with the origin shifting in k-space. In a ferromagnetic 
topological insulator thin film with out-of-plane magnetic 
anisotropy[I3I, the two surfaces are gapped, albeit with 
opposite mass. The mass domain wall created along the 
side walls of films then hosts a massless, one-dimensional, 
dissipation free edge state. If the chemical potential is 
tuned inside the magnetic gap of the surface states, trans¬ 
port is solely restricted to these ID edge states, and char¬ 
acterized by a quantized Hall resistivity p xy = h/e 2 and 
(ideally) vanishing longitudinal resistivity p xx = 0. This 
is known as the quantum anomalous Hall (QAH) state 
and was recently realized experimentally in thin films 
of Cr x (Bi,Sb) 2 _ x Te 3 [SI [B1 SJ at dilution fridge tempera¬ 
tures. 

Anisotropic magneto-resistance (AMR) PH E] is the 
response of the longitudinal resistance of a ferromagnet 
to the angle between the magnetization and direction of 
current flow. In this work, by tilting the magnetization of 
QAH insulator in-plane, we observe a novel, giant AMR 
effect that serves as a powerful quantitative probe of edge 
transport, even in temperature and chemical potential 
regimes away from perfect quantization (for instance, 
at high temperatures compared to the energy scale of 
the magnetic gap). Tilting the magnetization in-plane 
closes the magnetic gap and consequently destroys the 
edge states, leading to a crossover from dissipation-less to 


diffusive transport. We demonstrate that this crossover 
yields a sharp change in the longitudinal resistivity (p xx ) 
whose angular dependence is unlike classical AMR and 
show that modeling this dependence can be used to quan¬ 
tify edge state contributions to transport, and the extent 
of its inter-mixing with dissipative channels. 

The magnetic topological insulator devices used in this 
study were fabricated from 8 quintuple layer thick epitax¬ 
ial films of Cr x (Bi,Sb) 2 _ x Te 3 grown by molecular beam 
epitaxy on (111) SrTi0 3 (STO) substrates and measured 
in a vector magnet cryostat down to 280 mK (Figure 
la, see Methods for details). The large dielectric con¬ 
stant of the STO substrate at cryogenic temperatures 
allows effective electrical back gating for tuning the po¬ 
sition of the chemical potential[T8]. Three samples were 
studied with varying proximity to the QAH state (Ta¬ 
ble 1). In sample A, even at 280 mK, the AHE takes a 
maximum value p xy ~ 0.95h/e 2 (cr xy ^ 0.99h/e 2 ) and the 
longitudinal resitivity its minimum value p xx ~ 0.25h/e 2 , 
both at zero field, at charge neutrality (Vq 0 = — 100V). 
This non-vanishing p xx , accompanied by p xy taking near- 
quantized values, is indicative of dissipative channels co¬ 
existing with the ID chiral edge modes. p xx also displays 
a prominent hysteretic peak at magnetization reversal, 
and takes very large values approaching 260%. This 
large effect has been interpreted as enhanced backscat- 
tering between edge channels via domain walls at the co¬ 
ercive field, in a manner similar to the sharp p^ rise in 
the quantum Hall effect as scattering channels are pop¬ 
ulated (Figure lb,c). Below 3K near charge neutrality 
(Vg=Vg°), Pxx displays a metallic behavior as p xy ap¬ 
proaches quantization (Figure ld,e). This metallic be¬ 
havior is another signature of edge state transport and 
p xx recovers its insulating behavior when the chemical 
potential is tuned far out of the magnetic gap (see Sup¬ 
plementary information). Similarly, in the gate voltage 
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FIG. 1. (a) Cartoon of device and measurement scheme. 

Current flows in the x direction, and z is the direction per¬ 
pendicular to plane, (b) Near-quantized anomalous Hall effect 
at charge neutrality (Vg° = — 100V) and T = 280 mK, with 
p x y taking a maximum value of 0.95 h/e 2 at zero-field, (c) 
Corresponding longitudinal magneto-resistance, displaying a 
minimum p xx of 0.25 h/e 2 , also at zero-field, (d) and (e) 
Characteristic signatures of quantum anomalous Hall effect in 
the temperature and gate voltage dependence, respectively, of 
p xx (red) and p xy (blue) at zero-field (after magnetic train¬ 
ing). As edge-state transport dominates over diffusive trans¬ 
port, the decrease in p xx is complemented by a rise in p xy . 


TABLE I. Comparison of transport properties and AMR fit¬ 
ting parameters of 3 different Cr x (Bi,Sb) 2 - x Te 3 thin film 
devices with varying degrees of dissipation. 


Device p xy (/i/e 2 ) 

Pxx(0)(/i/e 2 ) 

V 

Rl 1 

AMR (%) 

A 

0.95 

0.25 

0.66 

0.62 

140.9 

B 

0.72 

0.55 

0.16 

0.54 

75.6 

C 

0.48 

2.79 

0.013 

0.025 

100.2 


dependence, as we tune the chemical potential towards 
the magnetic gap, p xx drops sharply, as p xy rises and 
finally flattens close to its maximum value. 

Having established the signatures of edge state trans¬ 
port, we now examine their evolution as the magnetiza¬ 
tion is tilted in-plane by an external IT radial field (Fig¬ 
ure 2). The temperature dependence of p xx in Figure 
2a shows a metal-insulator transition as the magnetiza¬ 
tion is tilted from out-of-plane to in-plane, irrespective 
of the azimuthal angle. Similarly, in the gate voltage de¬ 
pendence (Figure 2b), as the chemical potential is tuned 
into the magnetic gap, the drop in p xx is replaced by a 
sharp rise when the magnetization is in-plane. This is in¬ 
dicative of a crossover from dominant ID edge transport 
to diffusive transport, as the system transitions from a 
QAH regime to a gapless ferromagnetic topological in¬ 


sulator (when magnetized in-plane). Figures 2a and 2b 
already reveal a “giant” AMR effect when the chemical 
potential is in the magnetic gap, and this crossover is the 
origin of the novel angular dependence of the AMR that 
we now discuss. 

We systematically probe the crossover from edge- 
dominated to diffusive transport as we rotate a 1 T ex¬ 
ternal magnetic field from perpendicular-to-plane to in¬ 
plane. For a 1 T field, the AMR, defined as (p XXi max ~ 
Pxx,min)/Pxx,min takes a maximum value ~ 140%, and 
its angular dependence is completely inconsistent with 
the well-known cos 2 6 angular dependence seen in con¬ 
ventional ferromagnets. To explain the observations, we 
model the angular dependence using a simplified four- 
terminal Landauer-Buttiker formalism /j = T,jGij(Vi — 
Vj ) that assumes the coexistence of both dissipation-less 
edge channels and dissipative bulk and surface channels. 
The dissipative channels likely arise from the gapless, 
two-dimensional (2D) surface Dirac states of the side 
walls when the chemical potential is in the magnetic 
gappH- We note that these surface states are not gapped 
because the external magnetic field and the magnetiza¬ 
tion are always parallel to the plane of the side walls. We 
cannot exclude other possible origins of dissipation, such 
as the bulk bands, mid-gap impurity bands, and even the 
gapped surface states. These details are however not im¬ 
portant for the model. The conductance G tJ is given by 
Gij = Tij{h./e 2 ) where the transmission co-efficient Tjj 
takes the form I/,- = ijSij+i + tij. Here, r] is the con¬ 
tribution from edge modes that only transmit between 
leads j and j + 1, and t irj corresponds to dissipative con¬ 
tributions. From calculations explained in detail in the 
supplementary section, we derive the following expression 
for the magneto-resistance: 


eiTiiZ + Rdri + Rj 1 ’ 

where Rd is a measure of the “dissipative” resistance (and 
R^ 1 is the dissipative transmission coefficient). Physi¬ 
cally, the first term in the denominator is the conduc¬ 
tance of an inter-mixed channel, that constitues the edge 
states scattering into dissipative pathways. The second 
term though corresponds solely to the dissipative chan¬ 
nels. The transmission coefficient for the edge modes i] 
depends on the gap A induced by the magnetization. We 
expect ij to be suppressed from its maximum value of 1 as 
the magnetic gap is reduced A: this arises from enhanced 
inter-mixing with dissipative channels that accompanies 
an increase in penetration length A of the edge-state wave 
functions, where A = Hvf/A where vf is the Fermi ve¬ 
locity. We thus assume a phenomenological expression 
for the edge state transmission as: 

r? = 1 - e~ lo/x = 1 - e~ A/A ° (2) 

where lo is a length scale that quantifies the strength of 
the inter-mixing. This expression satisfies the physical 















3 



0 20 40 60 80 

V G - V G ° (V) 



0 50 100 150 200 250 300 350 

0 (degrees) 




FIG. 2. An external radial field of 1 T is used for all the data in this figure, (a) Transition from metallic to insulating behavior 
in the temperature dependence of p xx at Vg = Vg°, as the magnetization is tilted from out-of-plane to in-plane. The red, green 
and blue curves correspond to the x, y, and z directions respectively, (b) Gate voltage dependence of p xx at T = 280 mK when 
the film is magnetized along the x (red ),y (green) and z (blue) directions, (c) Out-of-plane anisotropic magneto-resistance of 
device A (light red, circles) at charge neutrality, as the field is rotated in the x — z plane. The angular dependence is unexplained 
by a conventional cos 2 9 fit (red, dotted line) and is understood by fits (red, bold line) that use an edge-dissipation mixing 
model. AMR data (circles) and corresponding fits (bold lines) from device B (light blue) and device C (light green), described 
in the text, (d) Cartoon of the device geometry highlighting field rotation from out-of-plane to in-plane, in the x—z plane, 
(e) In-plane anisotropic magneto-resistance of device A (light red, circles) at charge neutrality, as the field is rotated in the 
x — y plane. The angular dependence is explained by a typical cos 2 <j> fit (red, bold line), (f) Cartoon of the device geometry 
highlighting field rotation in the x — y plane. 


requirement at the two extrema: A —► 0 => 77 —>■ 0 
and A —► 00 ==> 77 —>■ 1. The angle dependence of the 
magneto-resistance then arises via the dependence of the 
magnetic gap on the perpendicular component of field, 
i.e. A = A m cos 6. A m is the size of the magnetic gap, 
when the film is magnetized out-of-plane. The physics 
behind the model is depicted in Figure 3. We now use 
Eq. 1 to fit our angular data with two free parameters: 
Rd and A' (defined as A' = A 0 /A m ). The model clearly 
fits the data very well, up to tilt angles ~ 70°, as seen in 
Figure 2b, confirming that the novel, giant AMR effect 
has its origins in the edge states. The model deviates at 
large tilt angles when the field approaches the in-plane 
orientation. In this regime, the dominant resistance is 
associated with the dissipative channels, whose angular 
magneto-resistance is not accounted for in our model (Rd 
is treated as a constant). 

Beyond explaining the angular dependence of this 
novel AMR, our model also presents a useful approach 
to quantify edge state contributions to transport in an 



FIG. 3. A pictorial representation of the physical origins of 
the giant AMR effect. The penetration length of edge states 
(depicted in dark blue) is dependent on the size of the mag¬ 
netic gap A, which in turn depends on the out-of-plane com¬ 
ponent of magnetization. As the external field is tilted in¬ 
plane, the size of the magnetic gap decreases (illustrated by 
the energy-momentum dispersions, in red). The penetration 
length of the edge states into the bulk (depicted in light blue) 
increases, as indicated in the middle panel, eventually lead¬ 
ing to complete mixing as the magnetization is tilted in-plane 
(last panel). 


edge-dissipation “mixed” QAH insulator. The fitting pa¬ 
rameter A' can be used to extract the edge transmission 
coefficient 77 = 1 — e -1 ' A . This is used to quantify and 
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compare edge transport contributions in two other sam¬ 
ples with varying degrees of dissipation, whose AMR is 
also shown in Figure 2c. The strength of the dissipa¬ 
tion is inferred by comparing the zero-field p xx and the 
magnitude of the AHE in these films. The angular depen¬ 
dence and fits to data from devices B and C also serves 
as a test of reproducibility for the giant AMR. Table 1 
summarizes the different device characteristics and the 
extracted fitting parameters, all at charge neutrality and 
at the same temperature. The extracted values of r] con¬ 
sistently characterize the degree to which the samples are 
close to the ideal QAH regime. Also, it is important to 
point out that the magnitude of the AMR effect consis¬ 
tently depends on the relative values of rj and R^ 1 (Eq. 
!)• 

As we mentioned earlier, conventional ferromagnets 
show an AMR described by a cos 2 <j> angular depen¬ 
dence. This behavior is indeed observed for in-plane 
measurements (Fig. 2e). Compared to the giant AMR, 
the in-plane AMR only has magnitude of ~ 10.7%, at 
charge neutrality (Figure 2e), and also shows a far weaker 
temperature dependence (difference between upper two 
curves in Fig. 2a). With the field and magnetization in¬ 
plane, the sample is expected to be in a diffusive regime, 
since tuning the surface states gapless destroys the edge 
states. More detailed mapping of the in-plane behavior 
is presented in the supplementary section. 

Finally, we present a more detailed picture of depen¬ 
dence of the giant AMR on the gate voltage and temper¬ 
ature. This is shown for device A in Figure 4. The effect 
is electrically tunable over an order of magnitude, and is 
largest in magnitude with the chemical potential in the 
magnetic gap, when edge state contributions to trans¬ 
port are maximum. As the chemical potential is tuned 
out of the magnetic gap, surface state and bulk dissi¬ 
pative channels are populated, reducing edge state con¬ 
tributions and consequently the magnitude of the AMR. 
This is captured in the gate dependence of the edge state 
transmission coefficient rj (Fig. 4b), extracted from our 
fits to the data of Fig. 4a. The qualitative behavior 
of R^ 1 is also reasonable, since the dissipative trans¬ 
mission coefficient increases as more electrons are accu¬ 
mulated. Similarly, the temperature dependence of the 
AMR (Fig. 4c) also captures the competition between 
the edge states and dissipative channels since our model 
quantifies edge transport in temperature regimes away 
from perfect quantization (Fig. 4d). While perfect Hall 
quantization has required dilution fridge temperatures, 
our results here demonstrate that the temperature scale 
for observing this novel AMR, with its origins in the edge- 
states, can be significantly higher. Furthermore, both the 
gate and temperature dependence of the fitting param¬ 
eters reveal the interplay between the edge states and 
dissipation channels: an increase in dissipative transmis¬ 
sion is accompanied by a reduction in edge transmission, 
and vice versa. 
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FIG. 4. (a) Out-of-plane AMR at T = 280 mK for a range of 
gate voltages. The AMR takes maximum values when the film 
is tuned to charge neutrality, (b) Gate voltage dependence of 
the edge-state transmission coefficient r/ (red) and the dissipa¬ 
tive transmission coefficient R^ 1 (blue) at 280 mK, extracted 
from fits to the data of (a), (c) Out-of-plane AMR at Vg° for 
a range of temperatures, (d) Temperature dependence of rj 
(red) and R^ 1 (blue) at Vg = Vg° extracted from fits to the 
data of (c). 


In summary, we have probed the evolution of edge 
states in a QAH insulator, as its magnetization is tilted 
in-plane, both by angular magneto-resistance measure¬ 
ments and by developing a model based on the Landauer- 
Biittiker formalism that accounts for co-existing dissipa¬ 
tive channels. We demonstrate a crossover from dom¬ 
inant ID edge transport to diffusive transport driven 
by the closing of the surface state magnetic gap. This 
crossover results in a novel, gate-tunable, giant AMR ef¬ 
fect that serves as a quantitative probe of the interplay 
between the chiral edge modes and dissipative channels, 
even in non-ideal QAH regimes. This methodology is rel¬ 
evant to ongoing efforts to raise the temperature scale of 
the quantum anomalous Hall effect [23] ■ Furthermore, as 
with recent advances that exploit the interplay between 
spin polarized states in topological insulators and ferro¬ 
magnetism PUH22] , the observation of a giant AMR effect 
may provide a useful basis for proof-of-concept devices for 
“topological spintronics.” 

Methods 

Material synthesis and characterization. The 
Cr x (Bi,Sb) 2 - x Te 3 thin films studied in this work were 
grown under ultra-high vacuum (low 10 -1 ° Torr) by 
molecular beam epitaxy on (111) STO . The source ma¬ 
terial was high purity Cr, Bi, Sb and Te thermally evap¬ 
orated from conventional Knudsen cells. For device C, 
e-beam evaporation of Cr was employed instead. Prior 
to growth, the STO substrates were subjected to an ex- 
situ high temperature (925 °C - 950 °C) anneal in an 
oxygen atmosphere to obtain a smoother substrate mor- 
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phology. The topological insulator film of device A was 
capped in-situ with a thin layer of room temperature 
deposited Al, which naturally oxidized upon exposure 
to atmosphere. Devices B and C were uncapped. The 
thickness of the topological insulator films was estimated 
by comparison with the growth rate of calibration sam¬ 
ples, whose thickness was measured using high resolution 
transmission electron microscopy and also confirmed by 
atomic force microscopy of the Hall bar step edge. 

Transport measurements All the devices used in this 
study were fabricated by mechanically scratching MBE- 
grown films into Hall bar geometries with typical lateral 
dimensions of around 1 mm. For example, device A has 
channel dimensions 1 mm x 0.5 mm. Since these Hall 
bars were mechanically scratched, the raw data was sym¬ 
metrized to account for geometric offsets. The angular 
measurements were performed in an Oxford Triton He 3 
vector magnet system with a radial field limit of 1 T and a 
base temperature of 280 mK. All transport measurements 
typically employed standard low frequency (19 Hz) lock- 
in techniques with source currents in the range 10-20 nA, 
and a Keithley 6430 source meter for back gating. 
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